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Abstract —Sparse code multiple access (SCMA) is a novel non- 
orthogonal multiple access technique, which fully exploits the 
shaping gain of multi-dimensional codewords. However, the lack 
of simplified multiuser detection algorithm prevents further im¬ 
plementation due to the inherently high computation complexity. 
In this paper, general SCMA detector algorithms based on Sum- 
product algorithm are elaborated. Then two improved algo¬ 
rithms are proposed, which simplify the detection structure and 
curtail exponent operations quantitatively in logarithm domain. 
Furthermore, to analyze these detection algorithms fairly, we 
derive theoretical expression of the average mutual information 
(AMI) of SCMA (SCMA-AMI), and employ a statistical method 
to calculate SCMA-AMI based specific detection algorithm. 
Simulation results show that the performance is almost as well as 
the based message passing algorithm in terms of both BER and 
AMI while the complexity is significantly decreased, compared 
to the traditional Max-Log approximation method. 

I. Introduction 

For the ability of supporting massive connections simul¬ 
taneously, Sparse code multiple access (SCMA) (TJ, a non- 
Orthogonal Multiple Access (NOMA) scheme, has been re¬ 
garded as a competitive candidate for Fifth Generation(5G) 
communication. Commonly, SCMA can be seen as a gener¬ 
alization of sparsely spread CDMA 0, with a few numbers 
of nonzero elements within a large signature length. When 
multiple layers are multiplexed with unique codebooks 0, the 
superposition scheme enable SCMA to benefit from shaping 
gain. However, to fully acquire these benefits, on the receiver 
side serious multiple address interference (MAI) is the main 
difficulty to process multiuser detection. The optimum maxi¬ 
mum a posterior (MAP) algorithm obviously shows the best 
detection performance with considerable calculation complex¬ 
ity. 

Meanwhile, some low complexity linear algorithm are pro¬ 
posed to handle this NP-complete problem a inevitably 
within some kind of performance loss. Especially, thanks 
to the sparse structure of SCMA spreading signature, the 
complex MAP formula 0 can be solved step by step within 
sum-product algorithm or message passing algorithm (MPA) 
0. Under an overloading system, detection complexity still 
increases with the number of users and size of codebooks. 
Thus, it is essential to take caution to simplify the algorithm 
and to find a tradeoff between complexity and performance. 

In this paper, inspired by the sum-product algorithm, we 
illustrate and prove how MPA is applied in multiuser SCMA 
detection and during the process of iteration, we then modify 
the message passing method to speed up the convergence in a 


reasonable way, which mainly make prior decisions on these 
high probability variable in advance. Furthermore, to reduce 
the calculation of product operation, we covert the whole 
algorithm into log domain and present the specific expression 
of symbols and bits soft message, which can be used in 
soft channel decoding conveniently. To eliminate the exponent 
operation, this work proposes two alternative methods based 
on the total least square criterion to approximate Jocabian 
logarithm 0 instead of simple Max-log 0. Moreover, in¬ 
spired by the concept of average mutual information in |9|, 
we take the SCMA scenario into account for further perfor¬ 
mance evaluation, proposing the average mutual information 
of SCMA (SCMA-AMI) and deriving the theoretical expres¬ 
sion of SCMA-AMI. Finally, the simulation results show the 
AMI gap between these simplified methods is negligible. The 
BER performance under both AWGN and Rayleigh channels 
are made a comparison between these algorithms, proposed 
methods present promising results superior to traditional Max- 
Log method. 

For the sake of clarity, the following notational convention 
unless noted otherwise. The sets of binary and complex num¬ 
bers are denoted by B and C. Upper-case calligraphic symbols 
denote sets, e.g., X. Scalars and vectors are represented 
by normal, bold fonts respectively. log(-)denotes the natural 
logarithm operation. 

II. System model 

We consider a general uplink SCMA system as shown in 
Fig(T| with J users spreading over K orthogonal resources 
(RE). In order to obtain higher spectrum efficiency, over¬ 
loading factor defined as /j = .// K is usually larger than 
one. On transmitting side, for each user j, j = 1,2...J, 
the binary information data stream {d}j = (dij,dk 0 j) 
are encoded by channel encoder with coding rate Rj then 
the corresponding output after code bit interleaver is {b}j = 
(bij, ..., b no j), where Rj = no/fco- These encoded bit streams 
are than mapped into AT-dimensional complex symbols, under 
the mapping relationship : / : B lo S 2 ( M ) Xj, Xj = 
f{b\j ,..., bo), Xj = (xij, ...XKj), where Xj belongs to C A 
with cardinality \Xj = M, D = log 2 M. Because of the sparse 
feature of SCMA, there are only Nj < K non-zero entries 
among a K dimensional constellation symbol x ? , the gener¬ 
ated codewords are then allocated on orthogonal resources, 
modulated by OFDM modulators. For simplification, here 
we assume that each layer corresponding to a specific user. 



Fig. 1. Illustration of a SCMA uplink system with J users spreading over K resources elements. 



Fig. 2. Factor graph representation of an SCMA system with J=6,K=4,N=2 


transmitting single bit stream. In fact, under prompt scheduling 
algorithm, more layers can be assigned to different users, while 
these subtle simplification has no essential influence on the 
subsequent analyzation. During a time slot, on receiver side, 
the received signal is the superposition of J users’ signal plus 
the ambient noise n, given by 

j 

y = ^diag( h J )x j + n , (1) 

i— 1 

where the received signal vector y = (yi,yx ) T , hj = 
(hij ...., hxj) T is the channel vector for user jth over K REs, 
and n ~ CJ\f( 0, NqT) is a K dimensional Complex Gaussian 
noise vector vector. 

The above uplink multiuser multiplex scheme can be de¬ 
picted in a bipartite graph G(V, £), which contains J variable 
nodes VN s, K function nodes FNs. Besides that, the edge 
between VNj and FNk means that resource k is occupied by 
user j the nodes. Let dj represents the set of FNs connected 
with the jth VN and dk has the similar meaning like dj. In 
an regular SCMA codebook, the degrees of VN s are equal 
to d„ = dj = N, while the degree of each FN is equal 
to dj = <9/,: which depends on the overloading factor. In 
Fig [2] shows an example of regular factor graph with J = 6, 
K = 4, N = dj = 2 and d v = 3. According to the above 
illustration and formula (1), it is easy to see the interference 
on fcth dimension as follows: 

Vk ^ ^ ^1'kjXkj + k 1 ,..., K. ( 2 ) 

jedk 

In next section, we will pay attention to the multiuser detection 
based on the elaborated system model. 


III. SCMA MULTIUSER DETECTION 

In this section, SCMA detector are elaborated in details 
based on message passing algorithm (MPA). Besides that, we 
improve and simplify the whole detection process. 

A. SCMA Detection with Sum-Product Algorithm 

The optimum detection criterion the object function is given 
by: 

K 

Xj = arg max V P({x 4 }) T] p(yk\{*i}J € dk), (3) 

a£ A zA A 
{xi}ex J fc=i 

x 3 =a 

where {xi} represent the set of constellation signals of related 
users. The key point is to effectively deal with the object 
marginalize product of functions. If the factor graph is cycle- 
free, the factor graph can be converted to an expression tree 
and each marginal posterior probability of VN can be obtained 
respectively. Further more, the redundant calculation of each 
VN can be shared simultaneously in a more effective method, 
which is proposed as Sum-Product algorithm or Message 
Passing Algorithm(MPA) @ in SCMA. 

Consider the bipartite graph depicted in FigU During each 
iteration, messages are first sent from VN s to FNs; Each FN 
then computes extrinsic messages and sends back to the VNs 
based on the previously received information. These VN-to- 
FN messages will then be used to calculate the new FN-to- 
VN messages in the next iteration. As \Xj \ = M, let Xj = 
1 > ; ■ • ■; ^.jM where Xj m — (•G jmi 

{ vj^ k (xj m )} represent the mth codeword message from VNj 
to FNk and {uj^^Xjm)} represents the message in the 
reverse direction at the fth iteration. These above messages all 
means the extrinsic message. Let dj/k denote the neighbor¬ 
hood of VNj excluding FNk and dj/k has similar meaning as 
above. H = {hi,..., hj,..., hj} represents the perfect channel 
estimation. The sum and product algorithm for computing the 
(approximate) posteriori distribution of all layers’ symbol are 
described as follows: in Algorithm 1. 

If the bipartite graph dose not exist any cycles, the iteration 
number T equals to the max degree of depth otherwise the 
correct object function can not be obtained, we can only attain 
an approximative result. 

In the final iteration, {Vj(xj m )} represent the approximate 
posterior probability of symbol x ;m , then each symbols of 

























































Algorithm 1 Sum-product algorithm for SCMA detection 
l: In put variable:y, H, X 

2 : Initial: P(x jm ) = 1/M, to = 1, = 1,J 

for j = 1, .J. k = 1 ,K and m=l,...,M do 
= 0 , - 1 
3: for t = 1 to T do 

firstly VNs —» PTVs 

for all j,k and m (if Edge(j,k) exists)do 

= P(Xjm) \xj m ). (4) 

l£dj\k 

To keep numerical stable, here normalize the data: 

M 

= VfH k (Xj m )/ ]T V&iXjm)- (5) 

m—1 

secondly TWs —>• CTVs 

for all />:, jand to (if Edge(j,k) exists) do 

| com | i 

Uk^j( x jrn) = ^-ea;p[- —(t/fc - h kj x k jm 

com ^ 

- E k «4r"V] n kE(x,„), (6) 

i£dk\j i€dk\j 

where com = ... x A) x ... x Xd f ~\ is the codewords 
”x” means Cartesian Product, wherei £ 9/c \j, \com\ = 
M d f~ l 

4: finnal result: for,/ = 1,J; m = 1 ,M 

Vj{ x jm) = P{xjm) U^j(xj m ). (7) 

kEdj 

5: Normalization 

M 

Vj (x^ m ) — kj' (Xjm,)/ ^ ^ Vj(p^jrri)i (^) 

m=l 


J layers can be determined at the same time. Here we also 
propose a simplified scheme: when observing the simulation 
process, we find that message passing algorithm is an adaptive 
algorithm, especially under cycled factor graph error detector 
probability will spread and influence other VN symbol’s 
judgement while some VN s’ symbol with higher confidence 
will quickly show the convergent tendency. Thus if in the 
foremost process of iteration, some VN show an obvious 
probability distribution, we can reasonable infer the correct 
transmit constellation on VN*. So we modify and simplify 
the above sum-product algorithm as follows: in Algorithm 2 

In this modified algorithm, if in the t-th iteration the symbol 
on some VN* is judged, when calculating {Uk^i}, where 
k £ dj*, VNi share the common FNk with VN*, the 
cardinality of set com decreases to M d f~ 2 , and {Uk^-j*} 
is no longer to be updated. Thus, the quantity of calculation 
exponent operations can sharply be curtailed. 


Algorithm 2 Modified Sum-product algorithm for SCMA 
detection_ 

l: main iteration : 

during the t th iteration 
for all m and j; 

V- emp {m) =P(x jm ) [J 

k£dj 

find m* = argmax[y^ emp (to)) in each layer j 
if max{y* emp (m*)) > a (a can be designed as 0.5) set 

( 1 TO = TO* 

Vj*(m) = Vj.-ykixj'm) = < _ , * ( 8 ) 

10 TO ^ TO 

wherefc £ dj* 

2 : other procedure is the same as sum-product algorithm 

when FNs connect the judged VNj*, com = ... x 
Xi x ... x Xd f -ithe cardinality of set com decreases to 
|com| = M^ df ~ 2 \ and there is no need to recalculate 
(xj* m ) any more because of the judgement of j* 
in advance. 


B. Simplified SCMA Dectection Based on LLR Type 

To manipulate the Sum-Product Algorithm practicably, the 
reduction of product calculation is essential. Thus, we utilize 
log domain to eliminate product operator. 

Based on the result in Section III-A: we soon get the 
following formula: let us fix a reference point x y i £ Xj, 
define {L^}. ,(x 7m )}, {L[*},-(x irn )} represent the log domain 
of {V} and {U} from formula (4) and (6) 
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L j^k 


(Xjm) 


V®, 


= log 3 
= log 


( x jrn ) 


(t) 

P ( X jm) 
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Cjm ~ jy (Vk hkjXkjm 
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final result 


Tj ( x jrn ) — log 


Vj( x jm) 

Vj( x j i) 


l P{ x jm) , T (T) i \ 77 

= l 0 g 744 y + E l€aAfc rn,j. 

( 12 ) 

Consequently, we proceed further and derive the output bit 
soft message LLR, as D = log 2 M bits are grouped together 















and then mapped to a constellation symbol, at the receiver side 
let A d represent the LLR of dth demapped bit. 

P{b d = 0|{V(x)}) 


A d =log 


P{b d = l|{V(x)}) 

E xe *C°) P ( bd = °I X ’ { V ( X )}) P ( X li y ( X )i) 

l0S E xP i-(i) p (.bd = 0|x,{y(x)})P(x|{y(x)}) 


= log 


E 




E 


(0)1/(x) E xG Af (0)e 


L(x) 


xg* 
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V(x) 
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xe* 


(i) e 


L(x) 


(13) 


where X^' denotes the constellation subset with the dth bit 
being b £ {0,1}. And then these bit soft message {A^j-can be 
further used in channel soft decoding. 

C. Low Complexity SCMA Detector 

As mentioned before, the difficulty to deal with the Log- 
MPA is in formula (6) while calculating exponent operation, 
we can regard the math problem as to process approximate 

\ com \ („) 

calculation: log c j™ and A,/ , where Jacobian logarithm 

,_, n —1 

[7] is feasible for this issue. 

log(e x + e v ) = f(x, y) = max(x , y) + log(l + e px ~ v ^) 

= max(x,y) + g(\x - y\). 

(14) 

Handling this equation we can calculate the above formula 
with iterative method 

N 

iog^V" = f{x n ,f(x n _i,f(x n _2,...,f(x 2 ,X 1 )...))). 

n= 1 

(15) 

We only need to approximate the second term of the right-hand 
part ofEq.d <?(•) in the iterative processing, because max(-) 
operator is easy to implement. Firstly completely eliminating 
the second part may be a easy choice, which is elaborated as 
’’Max-Log” g(z) = 0; 0 > 0 in H 

In this work, we use the total least squares criterion, 
choosing two curve types as prediction model to approximate 
the function g(-), the first non-linear curve is as follow: 

{ a 

- -he z < d 

b + z (16) 

0 z > d 

to optimize the total squared error between exact term and g i, 
we got optimal a = 1.0807, b = 1.1657, c = —0.1975, d = 5 
by means of numerical calculation and the mean squared error 
(MSE) is 5.2246 x 10" 5 . 

The second linear prediction model is given by: 


92{z) 


a(z — b) z < b 
0 z > b 


(17) 


In this model, the optimized parameter is a = —0.2614, 
b = 2.3555, the MSE between ,g 2 (-) and .9(') is 5.9254 x 10" 4 , 
obviously g \shows better fitting performance as shown in 
FigEl 



Fig. 3. The value of ln(l+exp(-z))and two approximate curves gl and g2. 


D. Average Mutual Information of SCMA (SCMA-AMI) 

Without loss of generality, we consider the Average Mutual 
Information (AMI) as a comparison criterion to evaluate the 
influence of approximation models a. 

In our sytem model, as bit streams are interleavered and 
then mapped to SCMA codebooks, on receiver side we can 
reasonably regard each bit level b d to a received y forms an 
independent channel. For a memoryless channel with perfect 
channel estimation, it can be assumed that input random 
variable X is equally taken from the points of the alphabet 
X, the conventional AMI is defined as IH: 

D 

= (18) 
d=l 

As we have stated in the previous chapter, SCMA is an over¬ 
loading multiplexing scheme, in such a case, it is reasonable 
to define SCMA-AMI as follows: 

D 

Iscma = At * Y\H ) =ti*J2 1{Bd - y ! p )’ ( 19 > 

d =1 


where [i = J/K, Bd E {0,1} is overloading factor, 

I{B d] Y\H) = E bd ^ h [I{b d] y\h)} 

= E bd ^ h [I(b d \h)-I(b d \y,h)\ 


— 1 Pb d ,y,h[log2 
= 1 ~ Eb d ,y,h[l°g2 
= 1 ~ Eb dt y t h[log2 


p(y,h) j 

p{b d ,y,h) 

E x g xP(v, h ,x) j 

Exe*«-d) p( b d, y, h, x) 
E x g;rP(y|M) , 
ExgA^PfelM) ’ 


( 20 ) 


Theoretically, if we use Maximum a Posteriori demapping 
criterion (MAP) or Message passing Algorithm (MPA) within 















cycle-free factor graph, we can derive the SCMA-AMI with 
formula (18) and (19). For further deduction, we define; 


t MAP 
^ d. 


log 


P(B d = 0\Y,H) 
P(B d =l\Y,H) 


since 


. 'E xeX mp(Y\x,H) 

( 21 ) 


Y l P(B d \Y,H) = l. (22) 


B d 


Following from (20) and (21), it is trivial to see: 

L MAP. . 

P(B d = b d \Y, H) = 6 “ \~ AP d) = P(B d \L™ AP ). 

1 + e L d 

(23) 


Using the definition of conditional mutual information and 
adding equation(21), (23) into the following formula, we can 
show that : 

HB d ; Y\H ) = I{B d ; L™ AP ) = H(B d ) - H(B d \L™ AP ) 

= 1 -E/ nb d ,L^ AP )\og nhd ^ AP) dLf AP : 

(24) 


As mentioned above, we got the theoretical method to obtain 
SCMA-AMI in formula (19) based on the output of optimal 
Log-MAP demapper. 

However, actually in SCMA system, the factor graph usually 
exists cycles, which results in inaccurate posterior output bit 
LLR X d m L^f AP and with different simplified algorithm 
quantity of estimation error will bring out different soft 
output. In detailed circumstance, we use I{B d \ X d ) to esimate 
the channel capacity and compare approximate SCMA-AMI 
between different demapper algorithm. 


I(B d -X d ) = h(X d )-h(B d \X d ), (25) 


h(X d ) = J p(Ad)log(-^-y)dAd, (26) 

J p(X d ,b d )log( ^ )dX d , (27) 

where/r(-) and h(- 1-) respectively denote the differential and 
conditional differential entropy. To obtain the PDF p(X d ) and 
p(X d \b d ), we can use the method of Mathematical Statistics 
to analyze the sample LLR X d , plot sample histogram to get 
approximate distribution data am. As demonstrated before, 
we insert formula (24) into formula (25) into (19) to get 
SCMA-AMI witin specific detector algorithm. In next section, 
we will give these simulation results. 

IV. Numerical results 

In this section, numerical results are presented to evaluate 
the performance of the proposed schemes. First, SCMA- 
AMI associated with different detection algorithms defined in 
Section III-D are shown in Figj4] for AWGN channel. This 
simulation is based on the basic parameters J = 6, K = 4, 


h(X d \B d ) = Y, 

B d 



Fig. 4. AMI Comparison of different SCMA detection schemes over AWGN 
channels 


N = 2, M = 4, p, = J/K{ connection relationship is shown in 
FigOJ- 6 layers SCMA codebooks are designed according to 
0- The SCMA-AMI loss engendered by simplified structure 
in algorithm 2 is negligible compared with original MPA in 
algorithm 1 as well as two approximation methods gl and 
g2. Meanwhile, the traditional Max-Log approximation shows 
obvious information loss especially in low SNR circumstance 
within almost 0.5dB deficiency. 

Second, BER simulation is carried out to evaluate the 
performance of these above algorithms. The regular systematic 
LDPC channel coding in China Mobile Multimedia Broadcast¬ 
ing (CMMB) systems is used in both AWGN and Rayleigh 
fading channel with the code rate R = 1/2 and 9216 coded 
bits per block. The channel coefficient h is set to a unit power, 
i.e., E(h 2 ) = 1. It is observed from Figj5] that in both of 
two types of channels the proposed three simplified method 
show almost the same BER performanc, the width of water 
fall regions are nearly 0.4 dB. Scarcely any losses in BER 
performance are gernerated by using simplified algorithm 1, 
gl and g2 approximate methods, while traditional Max-Log 
method DPC channel coding presents poor results, within as 
much as 0.5dB gap over AWGN channel and almost 0.25dB 
loss over Rayleigh fading channel. 

As what we have anticipated, three proposed methods 
present obvious advantages over traditional Max-Log methods 
with reasonable calculation complexity to some extent. 

V. Conclusions 

In this paper, we have proposed a simplified SCMA detec¬ 
tion algorithm via making prior decisions on some users with 
higher confidence, which reduces the complexity of further 
extrinsic information calculation on correlated edges. Further 
more, based on the logarithm expression derived in this paper, 
we approximated the iterative formula via curve fitting under 
the criterion of total least squares. On one hand, to analyze 
SCMA scheme in theory, the average mutual information of 
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Fig. 5. BER performance of SCMA detection algorithms over AWGN and 
Rayleigh Fading channels. 


SCMA (SCMA-AMI) has been proposed. Furthermore, the 
complexity of calculation is significantly reduced while the 
performance of proposed algorithm considerably outperforms 
that of using Max-Log approximation in numerical simulation 
results. In conclusion, proposed well-performed algorithms 
make an excellent tradeoff between the calculation complexity 
and performance in the SCMA system. 
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